Finite size effects in isobaric ratios 
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The properties of isobaric ratios, between nuclei produced in tire same reaction, are investigated 
using the canonical and grand-canonical statistical ensembles. Although the grand-canonical formu- 
lae furnish a means to correlate the ratios with the liquid drop parameters, finite size effects make 
it difficult to obtain their actual values from fitting nuclear collision data. 
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Central collisions between nuclei at energies starting at 
a few tens of MeV per nucleon lead the system to stages 
at which nuclear matter is hot and compressed [H-Q- If 
the coUision is violent enough, the subsequent expansion 
drives it to appreciably low densities and the nuclear dis- 
assembly takes place [H-Q- The properties of nuclear 
matter under these extreme conditions, i.e. its nuclear 
equation of state (EOS), play a very irnportant role in 
ruling the dynamics during these stages [j-Q- Thus, the 
study of the outcome of these reactions provides impor- 
tant insight into the EOS. 

Since the symmetry energy provide smaller contribu- 
tions to the equation of state, information about sym- 
metry energy is often obtained through ratios such as 
the isoscaling ratios constructed using the yields of nu- 
clei of mass and atomic numbers A and Z, respec- 
tively, produced in similar reactions, labeled '1' and '2', 
Y<^^^A,Z)/Y^^^A,Z) [3. Isoscaling has been found to 
be closely related to the difference between the chemi- 
cal potentials associated with the two systems. Thus it 
allows one to exploit its relationship with the symmetry 
energy [l^, [IH , although precise conclusions require very 
careful analyses [134l7| . 

Recently, it has been suggested by the authors of Rcf. 
[Tsj that important information on the liquid drop pa- 
rameters, including the difference in the chemical poten- 
tial of neutrons and protons, may be obtained by study- 
ing isobaric ratios between nuclei produced in the same 
reaction: 
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where I = N — Z . These authors offered an interpreta- 
tion based on a modified Fisher model [l^ to the different 
correlations studied in that work. More specifically, the 
yields Y{A, I) were calculated from this model and sev- 
eral manipulations, involving different ratios, were used 
to extract temperature dependent liquid drop model pa- 
rameters. The results were compared to the projectile 
fragmentation data reported in Ref. [20| . 



In this work, we resort to the standa rd g rand-canonical 
and canonical statistical approaches [2l| to investigate 
the extent to which reliable information may indeed be 
obtained from this analysis. Models based on these treat- 
ments have been extensively used in the last decades and 
have been very successful in describing many features of 
nuclear collisions [U [l^] . In the framework of the grand- 
canonical approach, the yields Y{A,I) reads: 
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where g stands for the spin degeneracy (taken as unit), 
\t ~ \/2-kW- IrrinT is the thermal wavelength, m„ de- 
notes the nucleon mass, /ip (/in) represents the proton 
(neutron) chemical potential, A/i = /i„ — [i^, fA,z{T) is 
the Helmholtz free energy associated with the fragment, 
and T is the breakup temperature. The free volume reads 
Vf = xVo, where Vq is the source's volume at normal 
density and we use x = 2 throughout this work. Upon 
inserting Eq. ^ into Eq. ([T]), one finds: 



lnRiI,I',A) = ^iI-I')- ^^'-'^^^-/-^-^'^^l (3) 

In order not to obscure the essential points of the 
present analysis, we adopt a simple prescription for the 
different contributions to the Helmholtz free energy: 



fAjiT) = ^Ba,i + fXAT) 



(i + x)i/3 



Z]/A, (4) 



with Zj = [A — /)/2. More specifically, the last term 
in the above equation corresponds to the Wigner-Seitz 
correction to the Coulomb energy [2^ and the binding 
energy 
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Ajj = 384 Z|,= 160 
Data: ^^Ca + ""se 
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FIG. 1: (Color online) isobaric ratios calculated for Odd A 
nuclei and for different system sizes in the framework of the 
Canonical ensemble at T = 6.0 MeV. The fit to the linear 
behavior shown in the figure gives Afi/T = 0.55, 0.95, 1.52 
and 2.00, while a*/T = 0.420, 0.280, 0.178 and 0.057 for 
the (48,20), (72,30), (96,40) and (384,160) nuclei, respec- 
tively. Results obtained with the grand-canonical ensemble 
(solid line) as well as ratios calculated with the experimental 
data of Ref. [20l | (full circles symbols) are also shown. 



is one of the simple formulae used in Ref. [24|, whose 
parameters read: a„ = 15.8 MeV, = 18.0 MeV, agym = 
23.5 MeV, and = 0.72 MeV. The internal Helmholtz 
free energy is the same employed in Refs. [l^, [13 
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for A> 5. It should be noted that this simple version of 
/* does not depend on /. The parameters, are the same 
used in Ref. i.e. Tc = 18.0 MeV, /3 = 18.0 MeV, and 
£0 = 16.0 MeV. Thus, Eq. ((3]) becomes: 



lni?(/,/',A) = ^(/-/')-^(Z|-Z|, 
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where we have introduced a* = ac[l — 1/(1 4- x)^/"^]. 

The equivalent formula derived in Ref. [1^ differs from 
this one by (small) terms involving the entropy of mix- 
ing. This small difference is slightly enhanced by factors 
associated with the pairing term of the binding energy, 
which we neglect in this work but is considered in Ref. 
p^ . However, since we confine our analysis to odd A, 
and focus on / = 1 and /' = —1, all these terms vanish 
and lni?(l,— 1,A) read: 
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where Z = {A ~ l)/2. The last term differs from that 
given by Eq. (8) of Ref. [l^ because those authors write 
the Coulomb energy proportional to Z[Z — 1) whereas 
we use Z^ . 

Figure [T] shows i?(l,— 1,A), calculated for Odd A 
nuclei, in both the grand-canonical and canonical ap- 
proaches, at T = 6.0 MeV. The grand-canonical yields 
are obtained by determining and /i„ from the con- 
straints Aa = Y.A,z^^A^z and Zq = Ea,z^^(^'^)' 
using Y{A,Z) given by Eq. Then, R{l,-l,A) is 

calculated directly from Eq. ([5]). The results are repre- 
sented in this figure by the solid line. By construction, 
itgives a* = 0.2208 MeV (the slope a*/T = 0.0368) 
12 [27|. Owing to the exponential relationship between the 
chemical potentials and the system size, A/x varies within 
1% in the considered mass range. Since the variation is 
very small, we use the ^^Ca value A/i = 12.5458 MeV 
(A/i/T = 2.091) in all the grand-canonical plots in this 
work. 

The results obtained with the canonical ensemble for 
different sources (Ao,Zo), i.e. (48,20), (72,30), (96,40), 
and (392, 160), arc also shown in Fig. [T] The linear de- 
pendence predicted by the grand-canonical formula, Eq. 
(|8|), gives way to curve lines which become more pro- 
nounced with smaller source size. The linear x depen- 
dence predicted by Eq. is a fairly accurate represen- 
tation of the actual behavior only for a; > 6, although it 
should also be valid for smaller x values. Furthermore, 
the fit parameters are strongly dependent on the system 
size and lead to values of A^/T and a*/T appreciably dif- 
ferent from the correct values. They seem to converge to 
the asymptotic values (predicted by the grand-canonical 
ensemble) only in the limit of very large systems. This 
' is illustrated by the symbols displayed in Fig. [51 which 

shows a*/T and A/i/T as a function of the source size 
Ao. 

This discrepancy is expected since, as discussed, for in- 
stance, in Refs. [2l|,[2^, finite size effects affect the yields 
predicted by the canonical ensemble so that it should be 
equivalent to the grand-canonical approach only in the 
limit when — >■ oo. Therefore, although the experi- 
mental -R(l, —1, A) ratios, also shown in Fig.[Tl exhibit an 
approximate linear behavior for x > 6, as asserted in Ref. 

(7) jl8| , the values of the parameters obtained from this anal- 
ysis are not clearly connected to the mass formula since 
the sources formed in actual experiments are not large 
enough to allow the finite size effects to be neglected. 
The projectile fragmentation data of ^^Ca -I- ^Be, which 
are depicted by the full circles in Fig.[Tl correspond to one 
of the systems studied in Ref. [l^. Note that the actual 
data do not exhibit a strictly linear trend as predicted by 
Eq. ([5]) and in Ref. [l^ . The staggering of the data points 
are not reproduced by any models and are very likely due 
to the deexcitation of the primordial hot fragments. Ig- 
noring this, the data follow the trend of the models with 
larger curvatures and therefore larger slope giving rise to 
slope parameter of a* /T = 0.75 and the offset parameter 

(8) of l^^/T = —1.7. From the present discussion, the ex- 



3 



^Be projectile fragmentation. 
- ^Be system as there are very 



tractcd parameters are expected to be strongly affected 
by finite size effects and, therefore, their physical inter- 
pretation is seriously compromised. It should be men- 
tioned that we observe the same qualitative behavior in 
the "^OCa + ^Be and ^^^Ni - 
We did not study the ^"^Ni 
few mirror nuclei produced in this case. It does not af- 
fect our discussion as the bottom line, i.e. the role played 
by finite size effects, does not depend on any particular 
source composition. Nevertheless, to avoid effects from 
different neutron to proton composition of the source, 
those used in the calculations in this work have the same 
N/Z ratios. 
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FIG. 2: (Color online) Size dependence of the fit parame- 
ters obtained with the canonical model. The dashed lines 
show the values predicted by the grand-canonical approach, 
i.e. al/T = 0.0368 and Afi/T = 2.091. For details, see the 
text. 



In conclusion, the linear behavior observed in the 
isobaric ratios studied in Ref. [H, for 2Z/A^^^ > 6, 
in the framework of a modified Fisher model can also 
be explained in terms of the standard canonical and 
grand-canonical statistical ensembles. However, the lat- 
ter strongly suggest that finite size effects prevent one 
from obtaining precise information on the nuclear prop- 
erties from such isobaric ratios. More specifically, the pa- 
rameters obtained in such analysis converge very slowly 
to the asymptotic values as a function of the system size. 
The effects are negligible only for system sizes which are 
much larger than those actually formed in the experi- 
ments. Thus, our results strongly suggest that it is not 
possible to safely ascribe a physical meaning to the pa- 
rameters obtained from such analysis. Furthermore, fi- 
nite size effects should not be neglected in the comparison 
of data to obscrvablcs obtained from the modified Fisher 
Model. 
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